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Abstract. The boundary double layer potential, or the Neumann-Poincare 
Oh operator, is studied on the Sobolev space of order 1/2 along the boundary, 

coinciding with the space of charges giving rise to double layer potentials with 
00. finite energy in the whole space. Poincare's program of studying the spec- 

trum of the boundary double layer potential is developed in complete gener- 
ality, on closed Lipschitz hypersurfaces in Euclidean space. Furthermore, the 
Neumann-Poincare operator is realized as a singular integral transform bear- 
ing similarities to the Beurling-Ahlfors transform in 2D. As an application, 
bounds for the spectrum of the Neumann-Poincare operator are derived from 
recent results in quasi-conformal mapping theory, in the case of planar curves 
with corners. 
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1. Introduction 

The boundary value problems of potential theory in K n are naturally cast as 
singular integral equations with unknown solutions supported by the boundary. 
When the boundary is singular, with finite vertex type corners as a mild example 
or it is a general Lipschitz hypersurface as an extreme example, notorious analytic 
complications arise. They were solved with ingenuity and perseverance during a 
century and a half of continuous struggle by many distinguished mathematicians 
and physicists. An early account of the fascinating history of this specific topic of 
potential theory is offered by the encyclopedia article by Lichtenstein |14) . 

Singular integrals related to layer potentials are thoroughly studied nowadays, 
first due to the central role they still hold and the source of inspiration they are 
in modern mathematical analysis, and second due to an array of new applications 
to continuum mechanics, field theory and engineering [2 r In this respect, the im- 
portance of the spectral analysis of the layer potential integral operators cannot 
be underestimated. The good news for the working mathematician is that much 
remains to be done on this front. For instance, the ubiquitous double layer poten- 
tial integral associated with the Laplacian, also known as the Neumann-Poincare 
operator, offers a very intriguing spectral picture. To be more specific, let T be 
a piecewise smooth planar Jordan curve with finitely many corners, endowed with 
arclength measure da and let K denote the Neumann-Poincare operator (see the 
preliminaries below for the precise definition); recent results of Irina Mitrea |16j 
imply a non-real spectrum of K on the Lebesgue spaces L p (T,da), 1 < p < oo, 
containing closed lemniscate domains, one for each vertex of T. On the other hand, 
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the spectrum of K on the Sobolev space H 1 ' 2 (T) 1 the natural home of charges 
of double layer potentials with finite energy, is real, and contained in the interval 
(—1,1]. A notable early contribution goes back to Carleman [5] who studied in 
great detail the resolvent of Neumann-Poincare's operator acting on the space of 
piecewise continuous bounded functions on V . Based on Poincare's original line of 
thought (see [H] for a modern account of it), one can safely state that considering 
the action of K on the space iJ 1 / 2 (P) stands aside as the most natural choice. 

Poincare's program, of estimating the spectrum of K : i? 1 / 2 (r) — > i/ 1 / 2 (r) 
via a Rayleigh quotient (balance) of outer and inner energies of the single layer 
fields generated by charges p G i? _1 / 2 (r), can be carried out in great generality. 
Specifically we will see its validity on closed Lipschitz hypersurfaces of R™, n > 2, 
in Section 3 below. Nonetheless, precise bounds for the spectrum of this operator 
are very scarce, see for instance [SJ. Two real dimensions are however special, due 
to the large group of conformal mappings. A groundbreaking observation due to 
Ahlfors [1 connects the spectrum of the Neumann-Poincare operator (acting on 
i? 1 / 2 (r)) to distortion estimates of quasi-conformal transforms, originally mapping 
the interior of T onto the exterior domain. A great deal of work has accumulated 
on the quasi-conformal estimates |12) . and as incomplete these results may be, 
they provide the most valuable source of bounds for the spectrum of the Neumann- 
Poincare operator in 2D. A second, and related, direction of research was advocated 
in 2D by M. Schiffer, who has obstinately returned during his career to the spatial 
interpretation of the Neumann-Poincare operator as a Beurling- Ahlfors transform 
acting on the Bergman space of the inner domain. Naturally, function theory of a 
complex variable plays a central role in Schiffer's context, see for instance [3]. 

The contents is the following. Section [2] develops Poincare's framework of inner 
and outer harmonic fields, in the case of Lipschitz hypersurfaces T in W l . The 
development is a matter of finding the proper spaces and formulations of results 
known for smoother varieties, and hence Section [2] is one of technical preliminaries. 

Section [3] generalizes results of |TT] to the Lipschitz boundary case. First by us- 
ing Plemelj's intertwining formula as the main ingredient in the similarity between 
the adjoint K* of the Neumann-Poincare operator and two different bounded, self- 
adjoint operators. These self-adjoint operators are in turn interpreted as angle 
operators between two canonical orthogonal decompositions of the space of har- 
monic fields of finite energy. In this way, Schiffer's and Poincare's ideas carry to 
any number of dimensions in the Lipschitz setting. A simple singular integral op- 
erator plays the role of the Beurling- Ahlfors transform, and Poincare's balance of 
energies is given a direct link to a self-adjoint angle operator. To give a flavor of 
the results, we state our theorem on the Beurling- Ahlfors transform in 2D. 

Theorem GEH Let Q C E 2 be an open and bounded Lipschitz domain with con- 
nected boundary and let Tq : L 2 (fl) — > L 2 (Q) denote the operator 



Then K* : H~ 1/2 {dn) -> 1/2 (dQ) is similar to Th : L 2 (ft) -> L 2 a {n), when the 
spaces are considered over the field of reals. Here Tnf(z) = To/(z). 

In Section 4 we exploit Theorem 13.61 and recent results in quasi-conformal map- 
ping theory to obtain bounds of the spectrum of the operator K in 2D, in the case 
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of domains with corners. An important point to be made is that SOT-type con- 
vergence arguments for the singular operator Tq, are exceedingly approachable, for 
example when approximating a domain Q with a sequence of domains Q n . To exem- 
plify, we give a simple and novel proof of Kuhnau's angle inequality; if a domain £1 
has a corner of angle Oj, then the spectral radius of Tq satisfies |er(Xh)| > \l — 8j/ir\, 
see Theorem 14.21 

In contrast with the spectral properties of K : L p — > L p [16], any non-essential 
point in the spectrum of K : H 1 / 2 — > H 1 / 2 is an isolated eigenvalue of finite 
multiplicity (see Corollary 13. 5|) . owing to the special symmetry features that K 
exhibits on H 1 / 2 . In modern computational applications involving K, or more 
general double layer potential operators such as those associated with the Helmholtz 
equation, such points in the discrete spectrum are easy to recognize. Furthermore, 
the resolvent can be controlled rather well in their vicinity, especially in comparison 
to points where there is a continuous contribution to the spectral picture. It is 
therefore of great interest to obtain information about the essential spectrum of 
the Neumann-Poincare operator. See for example [4], [9j. 

To state our main theorem, suppose that f2 is a curvilinear polygonal domain 
with interior angles smaller than ir such that its sum of exterior angles is less than 
2tt, with one angle regarded as negative. Then by the explicit construction of a 
conformal mapping from f2 onto a domain for which the spectral radius of Tq is 
known, we prove the following result on the essential spectrum of the boundary 
double layer potential. 

Theorem 14.61 Let Q be a C > a -smooth curvilinear polygon with Q < 6j < 7T for 
1 < j < N such that its angles satisfy 

N-l 

^2 0" - Oj) + 7T + 6 N < 2tt, 

possibly after a cyclic permutation of the vertex labels. Then 
\a ess (K)\ = \cr ess (T n )\ < max (1 - Oj/n) , 

l<i< AT 

where the Neumann-Poincare operator K is acting on 7? 1 / 2 (c3fi) . 

2. Preliminaries 

We begin by briefly introducing the Sobolev-Besov spaces and Lipschitz-Holder 
spaces that are necessary for our presentation. Let U C K™, n > 2, be an open 
Lipschitz domain with connected boundary. For a precise definition of a Lipschitz 
domain, see for example [2D]. ff 1 (t/) denotes the Hilbert space of functions V £ 
L 2 (U) such that 

H^lllfHc/) = H^||| 2(t/) + l|V^||| 2([/) < oo. 

H 1 (dU) may be instrinsically defined in a similar way, using the Lipschitz manifold 
structure of dU . For < s < 1, we obtain H s (dU) on the real interpolation scale 
between L 2 (dU) and H 1 (dU). Although it will not be used, we note that H s (dU) 
is a Besov space, 

\M 2 H s (au) ~ \\v\\ 2 L2(m) + f ^ X }~^ll da{ X )da{y), 

JdUxdU \ x y\ 
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where a denotes (n— l)-dimensional Hausdorff measure on dU. We define H~ s (dU), 
< s < 1, as the dual of H s (dU) under the (sesquilinear) L 2 -pairing, and by abuse 
of notation we always write 

9 if) = (f,g)L 2 {du)- 

Whenever 1 € H~ s (dU), H^(dU) denotes its annihilator, 

H s {dU) = {fe H S {8U) : (/, l) LHdu) = 0}, -l<s< 1. 

For d > an integer and < a < 1, we denote by C d,a (U) the space of d times 
continuously differentiable functions <f> : U — > C, such that d@(f> has a continuous 
extension to U satisfying an a- Holder condition, j3 = [fi\, . . . , j3 n ), j3j > 0, Yl Pj — 
d. That is, 

\d^(x)-d^(y)\<\x-yr, x,yeU. 

Next we recall a few conventions and facts pertaining to layer potentials. Let 
fi C M", n > 2, be a bounded Lipschitz domain with connected boundary. Let 
G(x, y) = G(x — y, 0) be the Newtonian kernel, normalized so that A x G(x, 0) = —5 
in the sense of distributions. Explicitly, 



G(x,y) = 



-w n 1 log|x-j/|, n = 2, 



where u> n is the measure of the unit sphere in R n . By the Neumann-Poincare 
operator, or the boundary double layer potential, K : _ff 1 / 2 (9f7) — > H 1 ^ 2 (dfl) we 
mean the operator 



Kf(x) = -2 f d ny G(x,y)f(y)da(y), 
Jon 



x G dtt, 

tan 



where n y denotes the outward normal derivative at y. Here the integral is intended 
in a principal value sense; its boundedness as an operator on L 2 (<9fi) was essentially 
proven by Coifman, Mcintosh and Meyer [8] , see also the classical work of Verchota 
|20) . The boundedness of K on i? 1 / 2 (9f2) follows by exploiting symmetry features 
of K and interpolation, see for example Lemma 13.11 By K* we always mean the 
adjoint with respect to the L 2 -pairing, so that, in our setting, K* is an operator 
acting on H~ l / 2 (dn). 
For x ^ <9f2 we write 



Df(x)= d ny G(x,y)f(y)da(y), x $ dCl, 
Jan 

and call D the double layer potential. For g 6 H^ 1 / 2 {dQ.) the single layer potential 
S is defined by 



Sg{x) = / G{x,y)g{y)da{y), x € K n . 
Jan 

Note that the operator S : L 2 (dfl) — > L 2 (dQ) is self-adjoint and bounded; 
furthermore we may assume that f2 is a normal domain in the sense that S : 
L 2 (dfl) — > H 1 (dfl) is a bijective bounded operator. By [20], this is equivalent 
to the existence of a function go e L 2 (dfl) such that Sgo\^ = 1. Note that 
ker(/ — K*) = <Cgo (see Lemma 13. ip . By a duality and interpolation argument, 
the extension S : H^ 1 / 2 (dVl) — > H 1 ^ 2 (dfl) is then also continuous and invertible. 
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If n > 3 every domain is normal, but when n = 2 there exist domains which are 
non-normal. However, if f2 is non-normal, then 

tQ = {tx : x £ fi} 

is normal (see [3D]) for every t > 0, t ^ 1. Since Kq and -fQn are unitarily equivalent 
and all our main results are invariant under such homotheties, we henceforth assume 
that f2 is normal. 

Denote the exterior of fl by fl e = £1 and by fj the space of harmonic functions 
h on U Q e with lim^^oo h(x) — and finite energy, 

\\h\\s3 = / \Vh\ 2 dx < oo. 

To ensure that fj is a Hilbert space we also require that if ft. 7^ and h e = h\a e = 0, 
then hi = ft|n is non-constant. 

Each element h £ Sj has an interior trace Txi n t h = Tr hi £ H 1 / 2 (dfl) and an 
exterior trace Tr ext h = Trh e £ H 1 / 2 ^®). By the classical Poincare inequality for 
bounded Lipschitz domains £/ and the fact that the trace Tr : iJ 1 (?7) — > H 1/>2 (dU) 
is continuous, we see that the interior and the exterior traces are continuous as 
maps from fj to H 1 ^ 2 (dfl). 

The trace normal derivatives d^h, <9° xt ft £ H~ 1 ^ 2 (dCl) are defined by duality 
and via Green's formula. To be more specific, 

(a**M>z»(an) = f(Vh,VVi)dx, (d°*%v) L * {dQ) = - f (Vh,VV e )dx, 

for u G iJ 1 /2(90) ) ^ g ffi(fi) with TrVi = v and K G H 1 ^) with TrV e = v. 
These definitions are meaningful. Suppose for example that TrVi — 0. Then there 
exist functions V n € C£°(Q) such that 14, — ► V< in [15] and a sequence of 

C°° domains f2j , compactly contained in and converging to il in an appropriate 
sense [20]. Then 




Letting n — > 00 we obtain that /q(V/i, VV<) g?x = 0, as desired. A similar argument 
shows that <9° xt ft is well-defined. On several occasions when we use Green's formula 
on Lipschitz domains in this paper, an argument like the one just presented is 
implicit. Note also that as maps, : fj -> H~ 1/2 {dVl) and <9£ xt : fj -> H- 1 / 2 {dQ) 
are continuous. 

One interprets an element h £ fj as a pair (hi,h e ) = (h\a,h\n e ), with the 
corresponding orthogonal decomposition fj = fji ® Sj e . We denote by Pi and P e 
the orthogonal projections onto fji and fj e respectively, so that Pi(hi,h e ) = (/ij,0) 
and P e {h % ,h e ) = (0, ft e ). 

Another natural orthogonal decomposition of fj holds, distinguishing among the 
single and double layer potentials. Specifically, let 

6 = {h E fj : Tr int h = Tr cxt h} 

denote the space of single layer potentials in fj, and let 

3) = {h £ Sj : B%*h = d^h} 

denote the space of double layer potentials. Then fj — 6©S) and we write P s and Pd 
for the corresponding projections. As expected we have that S : H~ 1 ' 2 (dO,) — > fj 
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is continuous and 6 = S(H- 1 / 2 (dCl)) if n > 3, while S : H 1/2 {dfl) -> £ is 
continuous and 6 = S(HQ 1 ^ 2 (dQ,)) when n = 2. Similarly, Z) : H^ 2 (dQ) — > Sj is 
continuous and S = D(H^ /2 (dn)) for all n > 2. Note that D(C) = C(1,0) is the 
subspace we subtracted from in order to have a norm. See Section 1 of for 
proofs which carry over verbatim to our Lipschitz setting for all the statements of 
this paragraph. 

We remark here that Poincare's inequality does not necessarily hold for the 
exterior domain, so that for an element h € fj, h e is not necessarily in H 1 ^,,). In 
spite of this, the exterior Dirichlet problem is well posed in the sense that Tr ext : 
$) e — s- H x / 2 {dVl) is continuous and invertible, where Sj e = P e $j. We have already 
seen that Tr ext is continuous, and its surjectivity follows by considering single layer 
potentials. It remains to check the injectivity. Suppose that Tr ext h e = 0. Then, 
for any V € fj e it follows that 

(d c n xt h e ,Tv cxt V) L 2 {an} = (Tr cxt h e , d^V) L ^ 9 n) = 0, 

and hence that also <9° xt ft. e = 0. But then (0, h e ) represents a function harmonic on 
all of E n and hence h e = 0. For if e C c °°(E n ), then 

f h e A0dx = - [ (Vh e , V4) dx = (dl xt h e ,^) L 2 ian) = 0, 

which implies that (0, h e ) is harmonic across dfl. This supplies the details for the 
terse argument presented in ([5], Proposition 4.1). 

For / e H^^idQ.) (/ _L 1 if n = 2) and g € H^ /2 (dQ), arguing with smooth 
functions and the continuity of operators involved, the well known jump formulae 
[2"0] for 5* and K take on the form 

Tr int Sf = Tr cxt Sf = Sf\ au , d^Sf = i(/ - K*f), 

dT'Sf = - K*f), Tr int Dg = l(-g - Kg), 

Tr cxt Dg = ±(g - Kg), d^Dg = d?Dg. 

3. The angle operators 

The present section is aimed at interpreting, up to similarities, the operator 
K* : H Q 1 ^ 2 (9J7) — > H Q 1 ^ 2 (dfl) as two different angle operators between the two 
orthogonal decompositions of the Hilbert space $). The first of these operators will 
be given a concrete realization as a generalized Beurling-Ahlfors singular integral 
transform, while the second will put Poincare's Rayleigh quotient of energies into 
its correct light in the case of a Lipschitz domain. The arguments essentially follow 
those of [11], with some additional technicalities, addressed in jS], arising from 
the fact that K is no longer a compact operator on the scale H s of Besov spaces, 
< s < 1. 

Lemma 3.1. For < s < 1 denote by W s C H s (dil) the subspace 

W s = {feH s : (f,g o ) L * = 0}. 

Then I — K : W s — > W s is invertible. In particular, for I — K : H s — > H s we have 
vaa(I-K) = W s . 
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Proof. By [20], I-K : W° -> W° and I-K* : L\ L% are invertible. Consider the 
single layer potential as an invertible operator S : L 2 — » H 1 . A simple computation 
shows that SL% = W l . From Plemelj's formula KS = SK* [6 we find 

I - K = S(I - K*)S-\ 

so that also I — K : W 1 — > W 1 is invertible. By a real interpolation argument it 
follows that I — K : W s — > W s is invertible. The last statement is now obvious 
since (I - K*)g = 0. 

□ 

Suppose that g e H 1 ^ 2 (dCl), or equivalently that Sg = — |(I — K)f for some 

1/2 

unique / 6 H (90), by the previous lemma and the fact that ker(7 — K) = C. 
By the jump formulae this equality means precisely that Df + Sg = in i7 e , which 
leads to the computation 



(1) 



\Df + Sg\\% = J (-±(f + Kf) + Sg^j d^(Df + Sg) da 

= - [ fW(Df + Sg) - d^(Df + Sg)) da = - [ fgda. 
Jdil Jon 



Summing up, J : H^ 1 ^ 2 —> H^/ 2 , J — 2(1 — K)~ 1 S is a strictly positive and 
bijective continuous operator. For future reference we also note that a similar 
computation shows that 

(2) (Df-Sg,Df + Sg) Si = - f Kfgda^(KJg,g) mm . 

As an operator considered on Lf,, J : L\ — » L\ is of course still strictly positive 
with dense range and hence has an injective square root \/~J with dense range. 



Lemma 3.2. The operator y/j extends to a bicontinuous bijection \fj : H { 



-1/2 



T2 



Proof. Noting that (g,f) H> (Jg,f)^ 2 is a scalar product for g e H 1 ^ 2 and 

1/2 

/ 6 H Q Cauchy-Schwarz inequality implies 

\(gJ) L A 2 = \{Jg,J- l f) L A 2 < (J9>g)L*{f,J- l f)v < \\VJ 9 \\Uf\\ 2 H1/ *, 

at least for g € Lq. Letting 1(f) = f gQ f da, this implies for / € H 1 / 2 that 

\{g,f)v\ = \{g,f-t(f)h>\ S W^Jgh'Wf - t{f)\\w,* < \\VJgh4f\\m/> 

We obtain for sufficiently smooth g, 

(3) \\g\\H-v><\\VJg\\zt. 

For h £ Lq one has the chain of inequalities 

\(Vjg,VJh) L 2\ < ||g|| H -l/2||J/l|| H l/2 < ||5(|| H -l/2||/l|| H -l/2 < ||g|| ff -l/2||\/j/l|| L 2. 

Due to the to the denseness of the range of vJ in h\ and a duality argument like 
above, one finds that yfj extends to a continuous operator \fj : H Q 1 ^ 2 — > Lq. We 
have seen in ([3]) that this operator is bounded from below, and since it also has 
dense range, the proof is complete. □ 



s 



KARL-MIKAEL PERFEKT AND MIHAI PUTINAR 



Consider for a moment the operator K\ : Lq —> Lq given by K\ — P h iK. Then, 
regarding the adjoint as an operator K\ : Lq —> L\ we infer K* = K*\ L 2 and 

K X J= JK*, 

the latter equation owing to Plemelj's symmetrization principle KS = SK* . Denot- 
ing by A : L\ — > L\ the operator A = \[jK*\[j , it follows that A is self-adjoint 
and that K* : H^ 1 ^ 2 — > H Q 1 ^ 2 is similar to A. 

— 1/2 

For g G H ,leth= —DJg + Sg G Sji as before, 
Sji = {h€Sj : h={h u 0)}. 
In the present language, (p} and ([2]) express that ||vJ<7Hz,2 = and 

((P d -P.)h,h)si = (AVjg,Vjg)^ o . 

Hence y/jg i— ¥ h is a unitary map of L\ onto Sji, giving that A is unitarily equivalent 
to the operator Pi(Pd — P s )Pi : iOi — i^i , since both operators are self-adjoint. The 
preceding computations and observations are summarized in the next proposition. 

Proposition 3.3. The operator K* (acting on H {dd)) is similar to the self- 
adjoint operator A = \fjK*\fj (acting on Lq(<951)). Furthermore, A is unitarily 
equivalent to the angle operator Pi (Pd — Ps) Pi acting on Sji — {h G Sj : h = (hi,0)}. 

Following the computations of in our Lipschitz setting we may now realize 
the angle operator Pi(Pd — P s )Pi as an operator acting on the M™-valued Bergman 
type space 25 (£1), 

Q5(0) = {Vu G L 2 (fl) : Au = 0}. 

Defining the operator 

IIn(Vu)(a;) = p. v. V x / V y G(x,y) ■ W y udy, x G 0, 
Jn 

we have the following result. 

Theorem 3.4. Let C K" be an open and bounded Lipschitz domain with con- 
nected boundary. The angle operator Pi(Pd~ P s ) Pi '■ Sji —¥ Sji is unitarily equivalent 
to the operator Bq : «8(fi) -> 93 (fi), 

B n =I + 2U n . 

In addition, the operator K* : H 1 ^ 2 (dfl) — > H Q 1 ^ 2 (dfl) is similar to Bq. 

Proof. Repeating the calculations of ([H], Lemma 5) verbatim, we have for / G 
ff 1/2 (dO) and g G H~ 1,2 {dVl) that 

p.v./ V v G(x,y) -V y Sg(y)dy = -Sg(x), x G Cl, 

and 

p.v. / VyG(x,y)-V y Df(y)dy = 0, xeft. 

Considering the unique decomposition h = Df + Sg of any h G Sji, the result follows 
via the unitary identification Sji 9 (hi, 0) i-> Vhi G B(Cl). □ 

Studying _ftT* acting on H 1 ^ 2 , rather than on H~ 1 / 2 , only eliminates the point 
1 from the spectrum of K* . 
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Corollary 3.5. The spectrum a{K) of K : H 1 / 2 — > H 1 / 2 coincides with that of 
Bq, except for the point 1. The essential spectra also coincide, 

a ess (K) = a ess (Bn). 

Furthermore, o~(K) C K and any point A £ a{K)\a ess {K) is an isolated eigenvalue 
of finite multiplicity, since Bq is self-adjoint. 

Remark 3.1. Chang and Lee [6] prove that cr(K) C (— 1, 1]. See also [9]. 

In two dimensions, n — 2, we may formulate Theorem 13.41 in terms of the usual 
Bergman space L 2 a (Vl) of analytic functions, 

L 2 a (Q) =L 2 (ft)nHol(D), 



its anti-analytic counterpart L 2 (tt), and the Beurling-Ahlfors transform. 

To make this formulation precise, we must consider S) and H 1 (dil) to consist 
only of real-valued functions. Then we may identify B(£l) with L 2 (fl), since any 
F € L 2 (fl) is of the form F = Whi = 2dhi for some h € fi. Note however that we 
are considering L 2 (fl) as a Hilbert space over the reals 

Let h — Df + Sg £ $)i for real- valued / £ H^ 2 and g £ H Q . Then a 
computation similar to that in Theorem 13.41 shows that 

P. v. - f d{ ?/ + 4 g) dA(0 = B(Df - Sg)(z), z £ il, 
* Jn (C - z) 2 

where dA denotes area measure. Note in particular that 

V f V,G(C !2 ) = i 7r ^. 

tt(C,~z) 2 

In conclusion, Theorem 13.41 can be restated as follows. 

Theorem 3.6. Let Q C R 2 be an open and bounded Lipschitz domain with con- 
nected boundary and let Tq : L 2 (fl) — > L 2 (Q) denote the operator 

T n /(z)=p.v.± / f^dA(C), f £ L 2 (f>), z £ Q. 

7T JO. (C - Z) 

Then K* : H Q 1/2 {dn) ->■ H~ 1/2 {dtt) is similar to Th : L 2 a {Vl) -> L 2 a {n), when the 
spaces are considered over the field of reals. Here Xh/(z) = Tq/(z). 

Note that the operator Tq is defined regardless of topological assumptions on fl 
such as boundedness and smoothness. It is also straightforward to check that if L 
is a fractional linear transformation, then Tq and T L ^ are unitarily equivalent. 
As in pTI] , we remark the following symmetries. Their proofs remain the same, 
except having to work with approximate eigenvalues rather than eigenvalues, noting 
that the spectrum of the symmetric operator Tq is equal to its set of approximate 
eigenvalues. 

Corollary 3.7. Let Q C R 2 be an open and bounded Lipschitz domain with con- 
nected boundary. Excepting the point 1, the spectrum o~(K) of K : H 1 ^ 2 (dfl) — > 
H 1 ^ 2 (dfl) is symmetric with respect to the origin. 

Corollary 3.8. Let C K 2 be an open and bounded Lipschitz domain with con- 
nected boundary. Then Q-{Tq) = a(TQ r ). 
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As the final matter of this section we shall consider an alternative angle operator, 
P s (P e — Pi)Ps ■ & — > &■ It seems geometrically plausible that this operator should 
have strong ties to the angle operator Pi(Pd — Ps)P%, and hence to the Neumann- 
Poincare operator. In fact, P s (P e — Pi)P s is intimately connected to Poincare's 
variational problem of inner and outer energies. To explain our point of view, let 
us briefly recall and combine some of the ideas of [9] and [IT]. The analogies with 
the previous work of this section should be clear. 

Let us restrict ourselves to n = 2 in our discussion. The higher dimensional 
cases are simpler. We will let y/S : H X ^ 2 — > W° play the role that yfj had 
earlier. Note that S : H~ 1/2 (dn) -> W 1 ' 2 is an invertible positive operator, so 
that y/S is invertible (cf. Lemma 13.21) . Define an operator B : W° — > W° by 
B = \fSK*\fS~ l . From Plemelj's formula KS = SK* we see that B is self- 

1/2 1/2 

adjoint, and by construction it is similar to K* : H — > H . Furthermore, 

— 1/2 

one easily verifies that for g € H Q (dfl) 

\\Sg\\% = {^Sg, ^Sg)L 2 (an), 

((P e - Pi)Sg, SgU = (BVSg, VSg) L 

We conclude that B is unitarily equivalent to P s (P e — Pi)P s : & — > ©. In particular, 
the spectral radius of K* may be computed via the numerical range of this operator, 



(4) \<r(K*\ H -y 2 ) 



\L 2 (n e ) w vu y\\L 2 (n) 

" iP i \^Sg\\l 2 ^ ) + \\VSg\\l 2m 
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which is exactly the supremum of Poincare's quotient. We state our conclusion in 
a theorem. 

Theorem 3.9. Let C M. n , n > 2, be an open and bounded Lipschitz domain with 
connected boundary. The operator K* : H-^ 2 (dfl) -> i?- 1 / 2 (9Q) (K* : H Q 1/2 -> 

— 1/2 

H if n = 2) is similar to the angle operator P s (P e — P%)Ps ■ & — > & acting on 
the space of single layer potentials. 



4. Spectral bounds for planar domains 

The points 

A g a(K) \{l} = a = o-(Th) 

are known as the Fredholm eigenvalues of J7, in the case that fi is a smooth bounded 
planar domain. The largest eigenvalue is often of interest. Theorem [376] allows us to 
define the largest Fredholm eigenvalue as |<r(Th)| for any simply connected domain 
f2 whose boundary is given by a closed Lipschitz curve in C* = CU{oo}. The most 
common definition of the largest Fredholm eigenvalue for a non-smooth domain is 
given by Poincare's extremal problem see [JJ]. Since they coincide for bounded 
Lipschitz domains by the results in Section [3] and both definitions are invariant 
under fractional linear transformations, Schober |19) . these two notions of a largest 
Fredholm eigenvalue coincide. We remark that the largest Fredholm eigenvalue is 
given by the numerical range of the symmetric operator To, 

(5) k(5h)l = ra= sup (2b/,/} LS(n) . 
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The largest Fredholm eigenvalue has been determined explicitly using quasicon- 
formal mapping techniques for certain domains, for example for regular polygons 
and rectangles sufficiently close to a square, see Kiihnau [T3] and Werner |21| . We 
shall make use of the following result of Krushkal [12] , proven by methods of holo- 
morphic motion. 

Theorem 4.1. |12j Let fl C R 2 be an unbounded convex domain with piecewise 
C l,a -smooth boundary, a > 0. Denote by < 9 < it the least interior angle made 
between the boundary arcs of dfl, taking into consideration also the angle made at 
oo. Then the largest Fredholm eigenvalue of £1 is 1 — 9/tt. 

Remark 4.1. Often the angle at infinity is considered to have the negative sign. 
For the purpose of this paper we will, however, consider it to be a positive angle. 
Despite this, it is admittedly instructive to think of it as a negative angle in some 
of the upcoming results. 

Denote by V the class of all domains ft described by the above theorem, as well 
as their fractional linear images L(f2). Hence, for every V G V we know that 

|ct(7V)| = max (1 - 9Jtt) , 

l<j<N J 

where 9i, . . . ,9n are the N interior angles of V. 

By a C 1 ' "-smooth curvilinear polygon f2 C R 2 , < a < 1, we shall mean an 
open, bounded and simply connected set whose boundary is curvilinear polygonal 
in the following sense: there are a finite number of counter-clockwise consecutive 
vertices (dj)f =1 C R 2 , 1 < N < oo, and C^-smooth arcs jj : [0,1] -> C with 
starting point <Zj and end point Oj+i such that dfl — Uj7j and jj and Jj+i meet 
at an interior angle < 9j+i < 2tt at the point <Zj+i. In this definition indices are 
to be understood modulo N. 

One major advantage of working with the 2D singular integral operator Tq is 
the easy accessibility to convergence arguments of SOT-type, for example when 
approximating a domain with better domains fl n . As an illustration, we gener- 
alize below, with a novel proof, the angle inequality of Kiihnau [13] , Theorem 14.21 
One should compare the approach below to the technical difficulties that appear 
when studying the convergence of the boundaries dfl n in order to prove that Kq^ 
converges to Kq SOT, explicitly carried out for instance in the excellent work of 
Verchota [20] . 

Theorem 4.2. Let fl C R 2 be a C l,a -smooth curvilinear polygon. Then 

|ff(5h)| > max 1 1 — . 

l<j<N J 

Proof. Without loss of generality we may assume that 

max II - 9i/n\ = |1 - 6 1 /n\ > 0, 

l<j<N 

a\ = and that 7i(0) > 0. It is then clear that when t — > oo, til converge as sets 
in C to the wedge 

W 6l = {re w :0<r<oo, 0<6>< 0{\. 

By Theorem O and Corollary ODH we know that \a(T Wei )\ = |1 - ^i/tt|, either W 9l 
or its complement being a convex domain with two corners of equal angles 9\ or 
2tt— 9\. For our technical needs, note that by mapping W$ 1 fractionally linearly into 
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the plane and applying the well known Markushcvich-Farrell theorem on polynomial 
approximation, one finds that any / £ L^Wg-J can be approximated by functions 
ft which are holomorphic on C* \ {zq} and square integrable in a neighborhood of 
oo, where C* = C U {00} and z e C is some fixed point outside Wg 1 . Note also 
that T t Q is unitarily equivalent to Tq for every t > 0. 

Denote by T : L 2 (C) — > L 2 (C) the Beurling transform, the unitary map of L 2 (C) 
given by 

Tf(z) = p. v. - [ ? i% dA(Q, f e L 2 (C), z e C. 
k Jc ((- z) 2 

For ft £ L^WgA, holomorphic on C* \ {z$} as described above, we have 

\\T t oh - T W0i h\\ L 2 {mnWgi) < \\T(x t n\w ei h )\\L 2 (c) + \\ T (Xw ei \tnh)\\L 2 (C) -> 

as t — > 00, since Xto\w 0l h — > and Xw 0l \tnh — > in L 2 (C). Here Xfl denotes the 
characteristic function of the set R. For similar reasons, 

In view of we may pick ft such that 

/ ftXV^ ft, cZA « \a(T Wei ) | . 

It then follows from the above computations that for t big enough, 

/ hT m hdAfa \a(T Wf>i )\, 

demonstrating that 

k(Th)| = kMI > \<r(TW~)\ = |1 - 0i/tt|, 
and hence completing the proof. □ 

We now turn to proving the bounds stated in Theorem 14.61 on the essential 
spectrum of Xh, and hence of K : H 1 / 2 — > H 1 ' 2 , for certain curvilinear polygons. 
The two lemmas below show that the essential spectrum of Xh is invariant under 
corner preserving conformal maps. 

Lemma 4.3. Let and Q>2 be C 1 ' -smooth curvilinear polygons with vertices 
( a j)jLi an d (bj)jLi respectively, such that the angle < Oj < ix made by dili at aj 
is equal to the angle made by at bj, for all! < j < N . Suppose that <f> : f2j — > Q 2 
is a biconformal map. It then automatically extends to a homeomorphism of£l\ onto 
0,2- Suppose that <fi( a j) — bj for all j . Then <ft G C 1,a (fl) and <fi'(z) ^ for z G Cl. 

Proof. By the classical Caratheodory theorem, <f> extends to a homeomorphism 
<fi : f2i — ► Sl2- By our assumptions, we may factor ip = ip2 ipi, where ^1 : fii — > D 
and tp2 '■ U — » ^2 are biconformal maps such that ipi(aj) = tj)^ (bj)- 

Fix for the moment a vertex aj of dfl± and let D e (tpi(aj)) be a sufficiently small 
closed disk of radius e centered at ipi(aj). For an appropriate branch of the power, 
it is easy to verify that the univalent map 

CH = (^r 1 H-%) 7r/9j 

on D then maps dH><~)D £ (ipi (aj)) onto a C^-smooth curve. By a local version of the 
Kellog-Warschawski theorem, see Pommerenke [18] , it follows that there exists an e 
such that C S C 1 ' a (Br\D s (ipi(aj))) and (' is zero-free close to ip\(aj). It follows that 
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C 1 also is C 1,Q -smooth in some relatively open neighborhood Uj = £(!£) n D e >(0) 
of 0, 

(6) C^W^i^ + a^^^). 
It follows that 

(7) M^ /7r + aj) = Ma 3 ) + zhj(z), 
where hj £ C°' a (Uj) and hj is zero- free. This leads to 

(8) i> 1 {z) = Tp 1 {a j ) + {z-a j yl (l !g j {z) 

holding in a neighborhood il± n D e »(aj), where gj is C 0, "-smooth, gj(dj) ^ 0. 
Similarly, by differentiating © we obtain in a relatively open neighborhood of a,j 
that 

(9) (z-a j ) 1 -^ d ^[(z) = e j (z), 

Cj being C 0, "-smooth and ej(a,j) ^ 0. Equations (jHJ and (|9]) express sharpened 
statements of (Pommerenke [17], Theorem 3.9). 

By arguing in the same manner for the map ?/>2 we find in a vicinity of ip\(a,j) 
that 

(10) (w-Ma j )) 1 - 9i/n ^M = f j (w), 

where fj is C°' Q -smooth and fj{4>i( a j)) ^ 0. 
In combining Q and (fT0|) we discover that 

tl/(z) = *e j {z)g j (z)W*- 1 fjW>i(*)), 
where A is some unimodular constant. We deduce that (f> is C ' "-smooth in a 
relatively open neighborhood of each vertex aj in Sli and 4>'{ a j) ^ 0. The validity 
of the same statement for each non- vertex point of dVli follows immediately by the 
local Kellog-Warschawski theorem. □ 

Lemma 4.4. Let Q,\ and VI2 be C la -smooth curvilinear polygons, and let <f> ■ ^1 — > 
0,2 be a biconformal map that extends to a homeomorphism tf> : fii — > Q2 such that 
4> G C 1 ' h (r2i) for some < b < 1 and <j)'(w) 7^ for w € Oi. Then Th 2 is unitarily 
equivalent to Tqj + K for a compact operator K : L^(Qi) — > L^(f2i), linear over 
the field of reals. 

Proof. Let the unitary map U : Lj^fla) — > L„(f2i) be given by 



Similarly, define t/ : L^(f22) — > ^a(^i) by Ug(w) — 4>'(w)g(<p(w)). Stokes' theorem 
and the change of variables C = fiiw) gives for a polynomial p that 



Tn 2 p(z) = lim — 
mo 2tt 



1 



2?r ./an 2 C - s 2?r Jaai </»(»?) 

Hence 



UT n2 p(w) = -- / ™WW drii 



27r iani </>(??) - 0(io) 
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and therefore the difference satisfies 



i 



(UTn 2 - T Ql U)p(w) = — I Up( V )K(n, w)dr h 

art! 



2tt 

where 



K(r],w) - 



<p(rj) — 4>(w) r] — w 
Since 4> g C 1,b (fli) and \<f>'\ > A > it is straightforward to verify that 

\K{v,w)\<\r]~w\ b -\ \d v K(r 1 ,w)\<\f 1 -w\ b - 2 , n,w£Q, 2 ,n^w. 
Hence, we may apply Stokes' theorem again to obtain 

(Z7To 2 - T ni U) P {w) = p.v. - [ Up{r j )d v K(r l ,w)dA(ri). 

That is, the difference UTq 2 — T^U is represented by a weakly singular kernel that 
is continuous off the diagonal, hence well known to be a compact operator (see e.g. 
[7]), completing the proof. □ 

Lemmas 14.31 and 14.41 combine in an obvious way. To prove Theorem 14.61 we 
require one more lemma, showing that the domains we consider can be mapped 
conformally onto domains in V in a way that preserves vertices. 

Lemma 4.5. Let il be a C 1 '" -smooth curvilinear polygon with < 6j < n for 
1 < 3 < N such that its angles satisfy 

N-l 

^2(^-0^+ it + 6 N <2vr, 
j=i 

possibly after a cyclic permutation of the vertex labels. Then there exists a curvi- 
linear polygon V G V and a biconformal map 4>:VL^V that extends to a C l > a {Q)- 
smooth homemorphism of Q, onto V, such that (j)'(z) ^ for z£!l. 

Proof. Let ipi : — > D be a biconformal map. By the classical Caratheodory 
theorem, it extends to a homeomorphism ipi : £1 — ► D. Let ip2 be the map defined 
on D by 

inspired by the Schwarz-Christoffel transformations. We will soon verify that 02 is 
univalent on ED with convex range. Based on this, one easily checks that dip2(P) C 
C* = CU {oo} is an unbounded piecewise real analytic Jordan curve having exactly 
N vertices, such that each vertex is given by ip2('>Pi(o,j)) for some j and the interior 
angle made there is equal to 8j. Note in particular that V ; 2('0i( a A r )) = 00 an d that 
2 (B) £ V. To obtain a bounded curvilinear polygon V G V, let V — -03 (02 
where 03 (z) = z -z ^ or a P om t z o 4- 02(B)- To finish the lemma, let <f> = 03°02°0i 



and apply Lemma 1473 
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To check the convexity of ip2, we mimic Johnston 10 and note that 

\ w-ipi{a N ) ^ w-ipi(aj) 

Since w h-> £ = maps ID biconformally onto 3?^ < 1/2 for z 6 c?B, it follows 
by our assumptions on the angles that 5ft ^1 + w ^j^j > for w € D. It is well 
known that this condition implies that is univalent with convex range. □ 

In combining the previous three lemmas we finally obtain our desired result. 

Theorem 4.6. Let Q be a C ,a -smooth curvilinear polygon with < 9j < tt for 
1 _ J _ A^ such that its angles satisfy 

N-l 

^2 0" - Oj) + n + On < 2tt, 
i=i 

possibly after a cyclic permutation of the vertex labels. Then 

\a ess (K)\ = \cr ess (T n )\ < max (1 - Oj/n) , 
i<j< n 

where the Neumann- Poincare operator K is acting on H 1 ^ 2 (dfl). 

As a finishing remark, we conjecture, partly based on numerical evidence, that 
I fess (?h) | only depends on the interior angles of a curvilinear polygon Q, and not 
on its shape. On the other hand, for the full spectrum, Werner [21] proves that 
|c(Tr)| = 1/2 for rectangles R sufficiently close to a square, but that \a(Tu)\ > 1/2 
when the ratio between the side lengths of R exceeds 2.76, suggesting the formation 
of eigenvalues outside the essential spectrum as a rectangle is elongated. 
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